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Abstract 



We have studied charmless hadronic decays of B mesons into two-body final 
states with kaons and pions and observe three new processes with the following 
branching fractions: B{B tt+tt") = (4.3l}-^±0.5) x 10"^, B{B K^ir^) = 
{U-Qtlitli) X 10-6, and B{B ^ K±^0) = (11.6l|0ti:^) x 10-^. We also 
update our previous measurements for the decays B — s- K^tt^ and B^ — > 
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CP violation in the Standard Model (SM) arises naturally from the complex phase in the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix Q. This picture is supported by nu- 
merous experimental constraints Q, as well as recent observation of direct CP violation in the kaon 
system but it remains an open experimental question whether this phase is the only source of 
CP violation in nature. Studies of the rare charmless decays of B mesons are likely to play an 
important role in constraining the CKM matrix and testing the SM picture of CP violation. 

Several approaches have been suggested to extract this phase information from measurements 
of rare B decays. Ratios of various B Kir branching fractions were shown |Q] to depend 
explicitly on 7 = Argiy*^^) with relatively modest model dependence. Within a factorization 
model, branching fractions of a large number of rare B decays can be parametrized by a small 
number of independent physical quantities, including 7, which can then be extracted through a 
global fit m to existing data. Finally, measurement of the time-dependent CP-violating asymmetry 
in the decay B^ vr+vr" can be used to determine the sum of 7 and the phase /3 = Arg{V^*^). In 
this case additional measurements of other isospin-related B ^ tttt processes are required to allow 
extraction of 7 -|- /3 |6| . 

In this Letter we present new measurements B ^ Kit and B — > vrvr branching fractions 
with significantly increased statistics, superseding results from our previous publication Q. In 
particular we present first observations of the long-awaited mode B tt~^it~ , as well as B ^ K^tt^ 
and B — > K^tt^ decays. 

The data used in this analysis were collected with the CLEO II detector at the Cornell Electron 
Storage Ring (CESR). It consists of 9.13 fb~^ taken at the T(4S), corresponding to 9.66M BB 
pairs, and 4.35 fb~^ taken below BB threshold, used for continuum background studies. 

CLEO II is a general purpose solenoidal magnet detector, described in detail elsewhere 
Cylindrical drift chambers in a 1.5T solenoidal magnetic field measure momentum and specific 
ionization {dE/dx) of charged particles. Photons are detected using a 7800-crystal CsI(Tl) elec- 
tromagnetic calorimeter. In the CLEO II. V detector configuration, the innermost chamber was 
replaced by a 3-layer, double-sided silicon vertex detector, and the gas in the main drift chamber 
was changed from an argon-ethane to a helium-propane mixture. As a result of these modifica- 
tions, the CLEO II. V portion of the data (2/3 of the total) has significantly improved particle 
identification and momentum resolution. 

Efficient track quality requirements are imposed on charged tracks, and pions and kaons are 
identified by dE/dx. The separation between kaons and pions for typical signal momenta p ~ 
2.6 GeV/c is 1.7 standard deviations {a) for CLEO II data and 2.0 a for CLEO II.V data. Candidate 

are selected from pairs of tracks forming well-measured displaced vertices with a tt+tt" invariant 
mass within 2a of the nominal mass. Pairs of photons with an invariant mass within 2. Sex of 
the nominal vr*^ mass are kinematically fitted with the mass constrained to the nominal vr*^ mass. 
Electrons are rejected based on dE/dx and the ratio of the track momentum to the associated 
shower energy in the Csl calorimeter; muons are rejected based on the penetration depth in the 
instrumented steel fiux return. 

The B decay candidate is identified via invariant mass and total energy of its decay products. 
We calculate a beam-constrained B mass M = \J E"^— p^, where is the B candidate momentum 
and E\^ is the beam energy. The resolution in M is dominated by the beam energy spread and 
ranges from 2.5 to 3.0 MeV, where the larger resolution corresponds to decay modes with a t:^ . 
We define /S.E = Ei + E2 — E]^, where E\ and E2 are the energies of the daughters of the B 
meson candidate. The resolution in AE is mode-dependent. For final states without 7r*^'s, the IS.E 
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resolution is 20 MeV (25 MeV in CLEO II). For modes with vr^'s the Ai? resolution is worse by 
approximately a factor of two and becomes slightly asymmetric because of energy loss out of the 
back of the Csl crystals. We accept events with M within 5.2-5.3 GeV and \I^E\ < 200 MeV. This 
fiducial region includes the signal region and a generous sideband for background normalization, tttt 
and Ktt signal events are distinguished both by dE/dx and IS.E observables. The Ai? distribution 
for B K^TT^ , calculated under the replacement of ttik by m,r, is centered at -42 MeV, giving a 
separation of 2.1a{1.7a in CLEO II) between B K^tt~ and B tt~^tt~ . 

We have studied backgrounds from b ^ c decays and other b ^ u and b ^ s decays and 
find that all are negligible for the analyses presented here. The main background arises from 
e+e~ qq (where q = u,d, s, c). Such events typically exhibit a two-jet structure and can produce 
high momentum back-to-back tracks in the fiducial region. To reduce contamination from these 
events, we calculate the angle 9s between the sphericity axis |Q of the candidate tracks and showers 
and the sphericity axis of the rest of the event. The distribution of cos 9s is strongly peaked at ±1 
for qq events and is nearly flat for BB events. We require | cos^sj < 0.8 which eliminates 83% of 
the background. Using a detailed GEANT-based Monte Carlo simulation |jl^ we determine overall 
detection efficiencies £" of 11 — 48%, as listed in Table |. Efficiencies include the branching fractions 
for Kg tt'^tt^ and vr'^ 77 where applicable. 

Additional discrimination between isotropic signal and rather jetty qq background is provided 
by the cosine of the angle between the candidate sphericity axis and beam axis (expected to be 
isotropic for signal, 1 + cos^9 distribution for qq background); the ratio of Fox- Wolfram moments 
H2/HQ |l^ (expected to be smaller for signal than for background); and the distribution of the 
energy from the rest of the event relative to the candidate's sphericity axis, as characterized by 
the energy in nine 10° angular bins. These 11 variables are combined by a Fisher discriminant 
technique as described in detail in Ref. |jl2| The Fisher discriminant is a linear combination of 
experimental observables J-" = X^i^i oaUi^ where the coefficients are chosen to maximize the 
separation between the simulated signal and background samples. 

We perform unbinned maximum-likelihood fits using AE, M, T , the angle between the B 
meson momentum and beam axis, and dE/dx (where applicable) as input information for each 
candidate event to determine the signal yields. Four different fits are performed, one for each 
topology [h^h^, /i^tt", h^K^, and K^tt^, referring to a charged kaon or pion). In each of 
these fits, the likelihood of the event is parametrized by the sum of probabilities for all relevant 
signal and background hypotheses, with relative weights determined by maximizing the likelihood 
function C. The probability of a particular hypothesis is calculated as a product of the probability 
density functions (PDFs) for each of the input variables. Further details about the likelihood fit 
can be found in Ref. ||l^. The parameters for the PDFs are determined from independent data 
and high-statistics Monte Carlo samples. We estimate a systematic error on the fitted yield by 
varying the PDFs used in the fit within their uncertainties. These uncertainties are dominated 
by the limited statistics in the independent data samples we used to determine the PDFs. The 
systematic errors on the measured branching fractions are obtained by adding this fit systematic 
in quadrature with the systematic error on the efficiency. 

Figure |l]a shows the results of the likelihood fit for B vr+vr" - and S ^ K^ir^. The curves 
represent the na contours, which correspond to the increase in — 21n£ by n^. Systematic uncer- 
tainties are not included in any contour plots. The statistical significance of a given signal yield is 
determined by repeating the fit with the signal yield fixed to be zero and recording the change in 
—2ln£.. We also compute from the PDFs the event-by-event probability to be signal or continuum 
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TABLE I. Summary of experimental results. The errors on branching fractions B are sta- 
tistical and systematic respectively. Reconstruction efficiency £ includes branching fractions of 
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background, as well as the probability to be /Cvr-like or vrTr-like. From these we form likelihood 
ratios, n,ig = {P:^ + PkJ/{P:^ + Pkn + + Pk. + Pkk) and = P^^/(P^\ + P^- Su- 
perscripts s and c denote signal and continuum background respectively. Figure ||b illustrates the 
distribution of events in TZsig (vertical axis) and TZt^ (horizontal axis). The cluster of events in the 
upper right corner is clear evidence for B tt+tt^. 

Figures ||(c-f) show distributions in M and AE for events after cuts on likelihood ratios TZsig 
and TZ-jt computed without M and AE, respectively. The likelihood fit projections for signal and 
background components, suitably scaled to account for the efficiencies of the additional cuts (50-70 
% for signal), are overlaid. Figure |2| shows the likelihood functions for the fits to — > K^-k^ and 
B /i=^7r°. Figure I shows M and AE distributions for B K'^tt^, B K^tt'^, and B K^-k^. 

We summarize all branching fractions and upper limits in Table |[ In addition to the first 
observations B —>■ tt'^tt~, B K~^7r^, and B — > K^tt^, we report improved measurements for the 
decays B — > K^ir^ and B — > K^tt^. The table also includes a range of theoretical predictions 
taken from recent literature p^. We see some indication for the decay B — > vr^vr'^ with the 



branching fraction of B{B — > tt^tt^) = (5.61^23 ^ ^■'^) ^ 10^^, but statistical significance of the 
signal yield is insufficient to claim an observation for this decay mode. We find no evidence for the 
decays B — > K^K~ and B —)■ K^K^, and calculate 90% confidence level (CL) upper limit yields 
by integrating the likelihood function 

where Cmaxi^) is the maximum £ at fixed to conservatively account for possible correlations 
among the free parameters in the fit. We then increase upper limit yields by their systematic errors 
and reduce detection efficiencies by their systematic errors to calculate branching fraction upper 
limits given in Table I. 

To evaluate how systematic uncertainties in the PDFs affect the statistical significance for 
modes where we report first observations, we repeated the fits for the h'^h^, h'^ir^ and K^tt^ 
modes with all PDFs changed simultaneously within their uncertainties to maximally reduce the 
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signal yield in the modes of interest. Under these extreme conditions, the significance of the first- 
observation modes 7r"'"7r~, K~^7r^ and K^tt^ becomes 3.2, 5.3 and 3.8 a respectively. We also 
evaluate the branching fractions with alternative analyses using tighter and looser cuts on the 
continuum-suppressing variable | cos^sj. These variations correspond to halving and doubling the 
background in the fitted sample. The changes in branching fractions under these variations are 
insignificant compared to the statistical error of our results. 

The ratio of the branching fractions B{B K^K^)/B{B tt^K^) can be used to estimate 
the size of final state interactions in charmless rare B decays [ 14 1 . Following the method outlined 
above we calculate B{B K^K^)/B{B t^^K^) < 0.3 at 90% CL. It has also been suggested |l§ 
to use the ratio of the branching fractions B{B — > K^7t^)/B{B tt^tt^) to estimate uncertainties 
in the measurement of the unitarity triangle parameter a = vr — /3 — 7 via B^{t) vr^vr^. We 
obtain B{B K^'k^)/B{B tt+tt') < 15 at 90% CL which implies that an error on a obtained 
from time-dependent asymmetry measurements of B^ vr+vr" can be as high as 60° [|l^ . 

In summary: we have made a first observation of B ^ vr+vr^; measured branching fractions for 
all four exclusive B Kit, including first observations of the decays B K^tt^ and B — s- K^ir^; 
obtained improved upper limits on i? — > vr^vr'^ and B — > KK modes. The hierarchy of branching 
fractions KK < vrvr < Ktt is obvious. 
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supported by the National Science Foundation, the U.S. Department of Energy, the Research 
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Foundation, the Swiss National Science Foundation, and Alexander von Humboldt Stiftung. 
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FIG. 1. Results for B K^tt^ and B tt+tt". Contours of the likelihood function versus 
Ktt and tttt event yield (a); likelihood ratios (b) — signal events cluster near the top of the figure, 

and separate into i^'-Tr-like events on the left and Trvr-like events on the right; beam constrained 
mass for i^vr-like events (c); AE for i^vr-like events (d); beam constrained mass for vr-Tr-like events 
(e); AE for TTTr-like events (f), with both tttt signal (dashed line) and Kir cross-feed (dot-dashed 
line) shown. 
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FIG. 3. Beam constrained mass and distributions for B — > K^ir^ (a-b), B K^tt^ (c-d), 

and B K^TT^ (e-f). 
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